Universal Features of Terahertz Absorption in Disordered Materials 
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Using an analytical theory, experimental terahertz time-domain spectroscopy data and numeri- 
cal evidence, we demonstrate that the frequency dependence of the absorption coupling coefficient 
between far-infrared photons and atomic vibrations in disordered materials has the universal func- 
tional form, C{lo) — a + Buj^ , where the material-specific constants A and B are related to the 
distributions of fluctuating charges obeying global and local charge neutrality, respectively. 



Since their name indicates a lack of order, it is perhaps 
unexpected that all disordered materials, regardless of 
their chemical composition, often behave very similarly in 
response to external probes. Low-temperature anomalies 
in thermodynamical properties of amorphous solids is one 
example [ij . Related to these is the pile- up of vibrational 
modes in excess of the Debye vibrational density of states 
(VDOS), g{Lo) ex w^, at typical frequencies ~ 30 cm~^, 
known as the boson peak (BP) (see e.g. 2] and Refs. 
therein). Interaction of photons with atomic vibrations 
in disordered materials is also expected to manifest uni- 
versal behaviourjespecially in the far-infrared (FIR) fre- 
quency domain [2|, 1 — 100 cni~^, where the atomic vi- 
brations are material-independent and resemble distorted 
sound waves. 

The simplest first-order perturbation process of pho- 
ton interaction with atomic vibrations is FIR absorption 
characterized by a linear absorption coefficient, a(w) = 
C{ijj)g{ijj), which is measured by IR spectroscopy — 
one of the most valuable experimental techniques to as- 
sess the vibrational properties of amorphous materials 
[J,l3 — and is found generally be a{uj) ex w^ in the FIR 
domain. The coefficient C(ijj) quantifies the degree of 
coupling between IR photons and atomic vibrations. It 
depends on the vibrational eigenmodes and the distribu- 
tion of atomic charges. The answer to the question about 
a possible universal frequency dependence of C{ijj) for all 
or most disordered materials in the FIR range has so far 
been uncertain, because: (i) C{uj) cannot be measured 
directly, and two independent experiments (for measur- 
ing a{uj) and g{i^)) are necessary; (ii) there is no rigorous 
analytical theory for a(tt') in disordered systems — previ- 
ous theoretical models concentrated mainly on the origin 
of the BP modes |a] within a particular soft-potential 



model or on the role of charge disorder in simple lattice 
models without positional disorder \j\. 

In this Letter, we demonstrate that, for disordered ma- 
terials in the FIR regime, C{ijj) indeed exhibits a univer- 
sal frequency dependence of the following form. 



C(w) ~A + Buj^ 



(1) 



where A and B are material-dependent constants. This 
law is valid in the low-frequency part of the FIR domain, 
namely for frequencies below the loffe-Regel crossover, 
ojiR, separating well propagating plane waves from 



w 
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strongly damped ones due to disorder-induced scattering 
[g|. In many materials, the loffe-Regel frequency is close 
to the BP 121 • We provide simple analytical arguments, 
confirmed by numerical analyses of molecular-dynamics 
(MD) models, and terahertz (THz) time-domain spec- 
troscopy (TDS) measurements of a(a;) for two glassy sys- 
tems, Si02 and AS2S3, and use independently available 
experimental VDOS data for g-Si02 [li| and g-As2S3 
[illE3 to extract C(w). 

Our analytical arguments are based on an analysis of 
the following general expression for the coefficient of ab- 
sorption of FIR photons by harmonic atomic vibrations 
obtained within the rigid- ion model jl!^ . 



C{u) ~ Co 



E 



ei{uj) 



(2) 



where Cq = li^'^njceJo , with mi and qi being the atomic 
masses and fixed atomic charges, 6^(0;) the component of 
the eigenvector of frequency cj corresponding to atom i, 
£00 stands for the high-frequency dielectric constant and 
n is the atomic concentration. A more general model 



with the fixed charges replaced by charge tensors [J, [ij] 
should be used for higher frequencies in order to describe 
properly the absorption peak positions and their relative 
intensities across the IR vibrational band, but, in the 
FIR range, the rigid-ion model is adequate, as we have 
checked for our density functional theory-based tight- 
binding (DFTB) [ll[l6J MD model of g-AsaSg |I3- 

In ordered systems, the eigenmodes are phonons and 
static atomic charges do not fluctuate and C{ui) is non- 
zero only for optic modes at the centre of the Brillouin 
zone. In disordered systems, structural disorder leads 
to a static charge transfer between atoms, i.e. to static 
disorder in the Qi, and to intrinsic disorder in the com- 
ponents of the vibrational eigenmodes. These two re- 
lated sources of disorder are encoded in Eq. (0) and are 
responsible for the peculiar behaviour of C{uj) in amor- 
phous systems. In order to demonstrate this, we use two 
known features about the structure of the vibrational 
eigenmodes in the FIR regime and about the distribu- 
tion of atomic charges in glasses. First, the disordered 
eigenmodes in the low-frequency regime can be approx- 
imately expanded in plane waves (see Ref. jg for more 
detail) characterized by a wavevector k and unit polariza- 
tion vector pk, i.e. ei{Lu) ~ J^u ^/rrn/rn ak(tf)pke'''^"-'""' 
with m = N^^ Tlii'^i ^^'^ ^i being the position vector 
for atom i {i ~ 1,...,N). The distribution of coeffi- 
cients ak(w) for low frequencies (at least below wir) has 
the shape of two relatively narrow peaks describing the 
disorder-induced hybridization between transverse and 
longitudinal acoustic phonons of the same frequency uj, 
with the peak positions scaling linearly with frequency 
[a [13. Second, it has been found for an ab-initio MD 
model of g-Si02 that the atomic-charge fluctuations pre- 
serve approximately local charge neutrality within Si04 
tetrahedra jj] . Similarly, we have found local charge neu- 
trality to be obeyed within AsSs pyramids in a DFTB 
MD model of g-AsaSa [13 ■ 

These two observations lead to Eq. (^). Assuming 
for simplicity that only a single plane wave with wave- 
vector k contributes to the disordered eigcnmodc (the 
presence of two peaks of finite widths in the distribu- 
tion ak(w) does not change qualitatively the results pre- 
sented below, at least for uj ^ wir), we can reduce the 
configurationally-averaged Eq. |j21l to the following form, 
(C(w)to/Co) ~ iV-i(|5|^), where S = J2^<l^e^"''■ Fur- 
ther analysis is based on the following rather general 
properties of the random sum S: (i) if the values of 
qi are not random, then the orientationally averaged S 
scales linearly with k for k ^ - the same holds for 
random but locally correlated values of qi (see below); 
and (ii) if the values of qi are random and uncorrelated, 
then S -^ const, for k ^ 0. It is convenient to split 
the atomic charges into two components, qi — qu + q2i, 
with gii({ri}) (depending on many atomic positions in 
a complicated fashion) representing uncorrelated charge 
components so that (quqij) = {qii){qij) = crfSij and e.g. 



{qiiS^^''^'} — {qii){e^^''^^} — 0. The random charges q2i 
satisfy local charge neutrality and can be imagined as re- 
sulting from charge transfers between nearest-neighbour 
atoms, i.e. q2i = J2i=ii ^Iji-^ where j runs over all nearest 
neighbours of atom i and Ag^i (= — Ag^j) is the charge 
transfer from originally neutral atoms j to i. In het- 
eropolar crystals, the values of Ag^i are finite and not 
random. In disordered systems, the values of Ag^i are 
distributed around mean values which do not necessarily 
coincide with those for their crystalline counterparts (see 
e.g. IJl). Such fiuctuations in Ag^i, and deviations of the 
means, are due to distortions in local structural units 
(e.g. the Si-O-Si bond angle in g-Si02) and the values 
of /S.qji are highly correlated on this length scale. The 
random charges qu are not correlated with q^j and thus 
{C{<jj)m/CQ) = iV-i(|S'ip + \S2?), where Sn coincides 
with S in which qi are replaced by qni- The first compo- 

nent, iV-^d^iH = A^"^ Ey ('Zi.9ij)(e*"'^> = <^l is in- 
dependent of k and thus of w. The second component can 
be evaluated in the bond representation, in which S2 = 
2i S (ij) ^9ije*'''^ sin (k-ry/2) with ry = Vj - r^ and 
Tij = {vj +ri)/2, and where the sum is taken over all the 
bonds {ij) in the system. In the FIR regime, k • r^ <C 1 
and hence 5*2 — fcX](ii) Aqye'''""'j(ik-rjj) = kS2- Conse- 
quently, the contribution to the coupling coefficient from 
correlated charges preserving local charge neutrality is 
N^^{\S2\'^) ^ fc^(|S'2p) oc co"^ (assuming linear disper- 



sion), where the function (|S'2| 



const. -I- 0(fc2) de- 



pends on precise structural details of the material but 
does not depend on k in the FIR range. Therefore, we 
have demonstrated that the frequency-independent part 
in Eq. (Q) is due to uncorrelated charge fluctuations while 
the quadratic frequency dependence results from corre- 
lated charge fiuctuations satisfying local charge neutral- 
ity. 

The range of validity of Eq. ^ is restricted at high fre- 
quencies by the requirement that the disordered modes 
should have acoustic-like character with linear pseudo- 
dispersion. Generally speaking, for lu > wjr, the disor- 
dered modes contain many plane waves with a wavevec- 
tor spread which is comparable with the mean value of 
k, and the dispersion law becomes modified, showing e.g. 
the onset of second sound for the longitudinal branch in 
g-Si02 la]. This means that the frequency dependence 
of C(w) can be modified in a non-universal manner for 
^ ^ ^iRj as we actually found for g-Si02 and g-As2S3. 
Dynamical polarization effects can also become signifi- 
cant for ui > wiR [19j. Eq. |^ is not restricted for low 
frequencies if possible low-frequency anharmonic effects 
are ignored Q- 

In order further to support the form of the univer- 
sal frequency dependence of the FIR coupling coefficient 
given by Eq. (^, we give the following evidence: (i) 
experimentally derived data for C{uj) in g-Si02 and g- 
AS2S3, and numerical calculations of C{uj) for (ii) clas- 



sical MD models of g-Si02 with fixed charges and for 
a— cristobahte (crystalline counterpart of g-Si02) with 
disordered charges and (iii) a DFTB MD model of g- 
AS2S3 and of its crystalline counterpart, orpiment, with 
disordered MuUiken charges. 

Experimental data for the FIR absorption coefficient 
of spectrosil WF g-Si02 and ultra-pure g-As2S3 at room 
temperature have been obtained by THz-TDS .2Qj in 
the frequency range 1 — 130 cm~^ using a TeraView 
TPI Spectra 1000 transmission spectrometer (experimen- 
tal details will be given elsewhere). Experimental data 
for a{uj) have been divided by the experimental VDOS 
obtained from inelastic neutron scattering [ifll, |ll| (see 
the insets in Fig. QJ and C{u!) thus obtained is shown 
in Fig. n Our data agree well with other experimen- 
tal FIR data obtained by THz-TDS or other techniques 
mm mil 13 and Ciuj) can be well fitted by Eq. (HJ 
in the frequency range of its applicability (see the dashed 
lines in Fig. ^. 
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FIG. 1: (Color online) Experimental TDS frequency depen- 
dence of the absorption coupling coefficient in the FIR range 
for (a) g-Si02 and (b) g-As2S3 (solid lines). Other experimen- 
tal data from Ref. |201 (o), Ref. (x), Ref. ^ {+) in (a) 
and Ref. [231 ("o'l in (b) are shown for comparison. The dot- 
dashed lines represent the numerical data obtained from MD 
models. The dashed lines show the fits of the experimental 
data by Eq.Elwith A = 4000 cm~^ B = 0.3 cm"^ for g-Si02 
and A = 1780 cm"^, B = 75 cm~^ for g-As2S3. The insets 
show the VDOS used to extract C{uj). The dashed line in the 
insets show the Debye VDOS, with the Debye temperature 
being 330 K for g-SiOa and 164 K for g-AsaSg. 

We have also calculated C{uj) for a classical MD model 
of g-SiO2 and a DFTB MD model of g-AsaSg [13 
(see the dot-dashed lines in Fig. QJ. The classical MD 
model of g-Si02 has fixed atomic charges qsi = 2.4 
and qo = —1.2, and thus C{uj) should not contain the 
frequency- independent component, (jS'ip), and should 
be proportional to w^ for lojI'Kc = a) ;$ cjir, ~ 35 cm~^ 
solely due to positional disorder. This region, unfortu- 
nately, is not available in our MD model due to prominent 
finite-size effects below 50 cm"-'^. However, in the range 



50 — 300 cm~^, we have found C(a;) oc u^ with /3 ~ 2 
which can be due to the existence of a peak-shaped spec- 
tral density and a linear pseudo-dispersion law in this fre- 
quency range |i]. The DFTB MD model of g-As2S3 incor- 
porates both charge and positional types of disorder and 
thus can show the crossover from a quadratic frequency 
dependence to the plateau behaviour with decreasing fre- 
quency (but this is also unavoidably masked by finite-size 
effects for a) < 20 cm^^) as seen in Fig. |2tb). Between 
the BP frequency, cDbp — 15 cm~^ 12] and 200 cm^^, we 
have found C(lo) cx oj^, with /3 ~ 1.7, which differs from 
that found for g-Si02 and might be due to non-linear 
pseudo-dispersion in this frequency range in g-As2S3 and 
possible dynamic-polarization effects [l9|. 
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FIG. 2: (Color online) Incoherent (a),(b) and coherent (c),(d) 
contributions to the coupling coefiicients for MD models of g- 
Si02 and a— cristobalite (a), (c) and of g-AsaSa and orpiment 
(b), (d) to the total coupling coefficient (dot-dashed fine). 
Solid lines represent crystalline models and dashed fines are 
for glassy ones. 

In order to unmask finite-size effects and reveal the 
role of charge fluctuations in the frequency dependence 
of C(a;), we argue that all the universal features of C{ijj) 
can be seen in lattice models of corresponding crys- 
talline counterparts incorporating correlated and uncor- 
related charge disorder. Many properties of disordered 
systems are rather similar to those of their crystalline 
counterparts (see e.g. [3, llg)- C{u}) shares such a 
similarity as well. In order to see this, it is conve- 
nient to split the expression for C(w) into incoherent 
(self-atom) and coherent (correlated-atom) components, 
C(w) = C''^™'^(w) + C™'^(w), where e.g. C'"™'>(cj) = 
CoX]j(9f/'^i)6i('^)- These quantities, calculated both 
for MD models of amorphous solids and their crystalline 
counterparts, show striking similarities (cf. the solid and 
dashed lines in Fig. |2a,b)), which suggests that the fre- 
quency dependence of C(lo) in amorphous systems can 
be mimicked by introducing correlated and uncorrelated 
charge disorder into crystalline systems (see also '7J). In- 



deed, we have found that uncorrelated disorder on its 
own, as expected, results in a frequency-independent cou- 
pUng constant (see the thin sohd Hues in Fig.EJ, the value 
of which depends on the variance of the distribution of 
qii, i.e. af. Correlated disorder on its own gives rise to 
the same frequency dependence (see the dashed lines in 
Fig. 0) of C{io) found above the BP in the MD models 
(cf. the dashed and solid lines in Fig.lHJ. The slope of the 
curves does not depend on the variance, a2, which only 
influences the intercept on the vertical axis. Incorpora- 
tion of both types of disorder (see the dot-dashed lines 
in Fig. O, with the plateau value coinciding with the 
parameter A found from the fit of Eq. (^ to experimen- 
tal data, reveals rather well the frequency dependence of 
C{uj) found in the MD models above the BP (cf. the 
dot-dashed and solid lines in Fig. O. Such a fit allows 
the widths of the distribution of uncorrelated charges, 
(Ti = ^yArn/Co, to be estimated, 171,3102 — 0.06 e and 
fi.AsaSa — 0.12 e; these values are in accord with those 
expected from the relative ionicities of the two mate- 
rials. Below 50 cm^^ (for g-Si02) and 20 cm^^ (for 
g-As2S3), the simulation data are not reliable due to 
finite-size effects. However, extrapolation to lower fre- 
quencies (larger sizes) is rather straightforward for lat- 
tice models with charge disorder. For example, in the 
case of g-Si02, the contribution from correlated charge 
fluctuations (the dashed line in Fig. Ola)) becomes in- 
signiflcant below the BP in comparison with the domi- 
nant frequency-independent contribution from uncorre- 
lated charge fluctuations (the thin solid line in Fig. ^a,) 
which should be extended to lower frequencies for larger 
models). Therefore, the sum of correlated and uncorre- 
lated charge fluctuations should result in a crossover to 
a plateau for frequencies below the BP, i.e. as in the 
experimentally observed behaviour of C{uj) (see Fig.^. 
In conclusion, we have presented an explanation for 
the universal frequency dependence of the coupling co- 
efficient for far-infrared photons with atomic vibrations. 
The coupling coefficient below the loffe-Regel crossover 
has two components. One is frequency independent 
and is due to uncorrelated static charge fluctuations 
caused by medium- and long-range structural irregulari- 
ties. This results in the quadratic frequency dependence 
of the absorption coefficient frequently observed in dis- 
ordered crystals and glasses. The other contribution de- 
pends quadratically on frequency and is caused by struc- 
tural disorder on the short-range (interatomic) scale, 
leading to static correlated charge fluctuations obeying 
local charge neutrality within structural units. For two 
glasses studied, g-Si02 and g-As2S3, we can conclude 
that, in g-Si02, uncorrelated charge fluctuations domi- 
nate through the whole frequency range below the loffe- 
Regel crossover (or boson peak) and result in a frequency- 
independent absorption coupling coefficient there. In 
contrast, in g-As2S3, uncorrelated charge fluctuations are 
less pronounced and correlated charge fluctuations pre- 
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FIG. 3: (Color online) Frequency dependence of C(a;) for MD 
models of a— cristobalite (a) and orpiment (b) with uncor- 
related (thin solid lines), correlated (dot-dashed lines) and 
both types of charge disorder (dashed lines). The uncor- 
related charges were drawn from normal distributions with 
standard deviations, asi/gsi = o"o/|qo| ~ 0.04 and ctas/^As = 
cs/l^sl = 0.06, and zero means. The correlated charges 
were taken from a normal distribution with mean values 
and standard deviations of: (g2i,Si) = 2.4, {q2i,o) = —1.2, 
{q2i,As) = 0.56, (q2i,s) = -0.37 and cr2,si/(g2i,Si) = 0.6, 
o"2,As/(g2i,As) = 0.6. The random values of q2i have been 
compensated by —q2i/Z placed on Z nearest neighbours in 
order to maintain local charge neutrality. The parameters 
of the distributions for uncorrelated charges were chosen to 
fit the experimental data for the plateau value (parameter A 
in Eq. (Q), and for correlated charges to fit the frequency 
dependence of C(aj) in the glassy MD models (solid lines). 



serving local charge neutrality become appreciable even 
below the BP and thus the absorption coupling coefficient 
exhibits an onset to an w^-dependence in this frequency 
range. 
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